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Abstract. We describe w-limit sets of completely positive (CP) maps over finite-dimensional spaces. 
In such sets and in its corresponding convex hulls, CP maps present isometric behavior and the states 
contained in it commute with each other. Motivated by these facts, we describe a quantization 
procedure based on CP maps which are induced by Markov (transfer) operators. Classical dynamics 
are described by an action over essentially bounded functions. A non-expansive linear map, which 
depends on a choice of a probability measure, is the centerpiece connecting phenomena over function 
and matrix spaces. 



1. Introduction 

In quantum mechanics there exists a great interest in problems in which the system and the 
surrounding environment are not isolated. Whenever we have that interactions between a system 
S and environment E cannot be ignored, we have an open system. The dynamics associated to the 
(reduced) system S is irreversible and incorporates noise and dissipative effects due to the presence 
of E. When considering open quantum systems, a starting point is the associated Master equation 
E21, 



(1-1) ^ = ~[H, P ] + ^[2L jP L* - {L*L 3 , P }], 

j 

where {A, B} = AB + BA, H is hermitian (coherent part of the dynamics) and Lj represents the 
coupling between S and E (Lindblad operators). Also, it is well-known that solutions for Master 
equations can be written in terms of completely positive maps (CP), $(p) = X^i^P^*' where the 
Vi are linear. In general, CP maps are not unitary and it is of interest to analyze the asymptotic 
evolution of the associated dynamics. This has been done for finite-dimensional CP maps in |13j . 
in the case of mixed-unitary operators, and [9J for bistochastic operators. A related question on 
CP dynamics is to describe the cu-limit sets, that is, the set of density matrices which are limit of 
$ ni (po) for a fixed density po, and for some subsequence {n^}. There are well-known results on limit 
sets coming from the theory of non-expansive maps over Banach spaces, so it is natural to pursue 
such description on CP maps, since these are non-expansive with respect to the trace distance, 
d(p, a) = \tr\p — a\, see [13] . 



Motivated by certain limit results, we consider a classical map which acts on densities (i.e., positive 
functions on a compact space). Assume such dynamics is described by a linear, positive map. We 
ask the following question: is there a CP evolution on some matrix space which can be seen as the 
quantum counterpart of the classical evolution? In this work we describe a quantization procedure 
based on CP maps. We will require that fixed points of the classical dynamics are taken to fixed 
points in the quantum setting, and similarly to w-limit sets. The procedure connects classical and 
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quantum dynamics in the following way: 
(1.2) 



L°°{E M -p) — ^ L°°{E M -p) — ^ 



I n ®M n > /„®M n y ■■■ 

The set L°°{E M ] p) denotes the set of essentially bounded functions with respect the measure p over 
the space Em of state functionals over M = M n (C), the order n matrices. The upper line represents 
the classical evolution and the operator L is the transfer operator. Its dual, acting on probability 
measures, is the associated Markov operator. The lower part contains A, which is a CP map to be 
determined. The map k m is defined as follows. Let p G M p (Em), a probability measure over Em 
with barycenter p. Then / G L°°{p) (->■ ^ -^i ® -^fn(C) is defined by: 

(1.3) (O p ,^(/)tt,(A)O p ) = / f(u')u)'(A) dp(u') 



Such map is linear and non-expansive [3] • We choose a faithful state p G M n and set the GNS triplet 
as (Hp, n p , Hp) = (C n <g> C n ,M„ (g) /„, |a/p)). By this we mean that ir p (M n ) = M n <g> J n , where if 
P = Y.i r i\ r i)( r i\ e M n, then define 

(1.4) \VP) = J2 v^l r *> ® l r *> e £H ® C " 

i 

For a certain class of measures (orthogonal), Tomita's theorem [3] states that is a ^-isomorphism 
into I n ® M n . But we note that the description presented in this work is also of interest even when 
the measures involved are not orthogonal. Given L a transfer operator, a positive function / and 
a probability measure /i, we wish to find A = A(/x, /) such that the above diagram commutes. We 
note that Em = Em u , the space of state functionals over C n , can be identified with the space D(C n ) 
of density matrices in M n (C). We prove in section [3J 

Theorem Q Let L : C{E Mn ) C(E Mn ), L{g){p) = £LiP^(p)), k < n\ Pl > 0. 
(1) For a given / > and for every p = Ylim ^m$u> m , A m > 0, oj, with support on the diagonal 
units of M n (C), there exists a solution for A o = Kp(Lf), given by A(p) = Yli=i ViP^C ■> 

Vi = diag(t>i, ...,<), and 



;i.s) ki = ^ 



Pit Em^/i^Wl^f^ii) 



fc, j = l,...,n 



Em ^mfi^mj^miEjj) 

(2) If L/ = Ylii"Pif{Fi) with Fi = Ui ■ U*, Ui unitary, there exists Vi which satisfies A o k m (/) = 
K/j,(Lf) for every / G L 00 if and only if the cjj are all fixed points for $(p) = YliPi-^iip)- 

Among the known connections between Markov operators and quantum systems, we recall that a 
density matrix p is a fixed point for a mixed unitary channel A if and only if p is the barycenter 
of a measure which is invariant for the associated Markov operator [TU], [TH]. We are interested in 
further developing such connection, having in mind applications to quantum information theory. In 
this article we focus on finite-dimensional quantum channels. 

The structure of this work is the following. In section 12.11 we review basic facts on completely 
positive maps. In section 12.21 we review w-limit sets, results on non-expansive maps, and describe 
the related results for CP maps which are of our interest. Several of the results in that section are 
easy adaptations of well-known theorems on non-expansive maps on Banach spaces and the proofs 
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are briefly described in the appendix. In section 13.11 we review some facts on measure theory and 
representations. In section 13.21 we review the transfer operator, which represents the action on the 
space of functions. In section 13.31 we describe the quantization procedure and conditions for its 
existence. We conclude in section [4] with some open questions on a classical commutative diagram 
and its relations to the construction made in this work. 

2. Evolution of classical and quantum dynamics 

2.1. Completely positive maps. We make a brief review on completely positive maps [14]. For 
simplicity we will often restrict ourselves to the space M n (C) of order n matrices. Let Tl be a Hilbert 
space, let B(H) be the set of bounded linear operators over 7-L. If dimTi < oo we will take H = C n , 
and in that case B(C n ) = M n (C). A matrix A : C n ->• C n is positive, denoted by A > 0, if 
(Av, v) > 0, for all v G C n . This is also equivalent to any of the following facts: a) all its eigenvalues 
are nonnegative; b) there exists a matrix B such that A = BB* . In particular, every positive matrix 
is hermitian. We say p G M n (C) is a density matrix if p > and tr(p) = 1. 

Let $ : M n (C) — > M n (C) be linear. We say $ is a positive operator whenever A > implies 
<&(A) > 0. In physics problems, we are mostly interested in <3> that are not only positive on M n 
(system A), but also that I £g> $ is also positive, where / is the identity operator on any other choice 
of component (system B). More precisely, define for each k > 1, : M k (M n (C)) — > Mj.(M n (C)), 

(2.1) a> k (A) = [$(Ai)L A e M k (M n (C)), Aij e M n (C) 

We say $ is A;-positive if $fc is positive, and we say $ is completely positive (CP) if is positive 
for every k = 1, 2, 3, ... . It is well-known that CP maps can be written in the Kraus form [14] : 

(2.2) $( p ) = J2V l pV: 

i 

Conversely, any such operator is CP. An operator of the form 

k 

(2.3) ${p) = Y lPi U i pU;, pi>0,J> = l 

i=l i 

where the U{ are unitary is called a mixed-unitary operator or a random unitary operator. 

Such operators are unital ($(/) = /) and trace-preserving. It is well-known that p is fixed for a 
mixed- unitary operator if and only if it is fixed for each Fi(p) = UipU* [SI [32]. 

2.2. w-limit sets. Let D(C N ) denote the compact, convex space of density matrices over C . We 
identify it with the space of state functionals over C . We consider a CPT (trace-preserving) map 
$ acting on D(C ). Define, for p G D(C ), the orbit through p, 

(2-4) 7 (p) = U n > 1 $ ri ( / 9), 

and the fi-limit set of p: 

(2.5) Q(p) = {r] G D(C n ) : r] = lim ^(p),^ -> oo as i -> oo}, 

i— >-oo 

Define the attractor of $ to be £7$ = U p fi(p). Other definitions are the following. 
Definition. 

(1) We say that a subset B C D(C N ) is positive invariant under $ if $™(_B) C 5 for all n > 1. 

(2) We say that B is minimal under $ if for any r\ G B = 7(77), that is, if B is the closure of 
the orbit of any of its elements. 
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(3) We say that B is strongly invariant under $ if for every n > 1 we have that $ n is a 
homeomorphism of B onto B. 

We describe some results which allows us to compare convex sets in the classical setting with 
the corresponding objects after the quantization procedure we will perform. Most steps are due to 
results proved in the setting of non-expansive dynamical systems. The main ideas are well-known 
[U E], but for completeness we describe the proof and simple adaptations needed for our setting 
in the appendix. Our analysis is mostly relevant for CP maps acting on several qubits where, for 
instance, we obtain periodic asymptotic dynamics, such as the two-qubit CNOT channel [T3]. When 
the asymptotic limit is trivial (e.g. amplitude damping channel for one qubit, which has the north 
pole of the Bloch sphere as limit for all density matrices), many of the results presented are trivial. 

Proposition 2.1. For all p G D(C N ) we have the following, a) fi(p) is minimal under $. b) For 
every k > 1, $ fc is an isometry on £l(p). c) fi(p) is strongly invariant under $. 

Recall that the convex hull of a set S, denoted co S, consists of all the convex combinations of 
points from that set. Also recall that by Caratheodory's theorem, we have that in d- dimensional 
Euclidean space, every point in co S is a convex combination of d + 1 points from S. The convex hull 
of a finite set of points is called a polytope. We will be interested in the convex hull for omega-limit 
sets associated to CP maps acting on finite-dimensional spaces. 

Lemma 2.2. If the restriction of $ to a subset B of D(<C N ) is an isometry then the restriction o/$ 
to co B is also an isometry. 

Proposition 2.3. a) The set co Q(p) is strongly invariant under $. Also there is an affine group 
of isometries on the closed linear variety spanned by Q(p) which coincides with $ on coQ(p). b) 
co fl(p) contains exactly one fixed point po o/$, given by 



Proposition 2.4. a) There exists a non-expansive projection r onto Q$. b) We have that Q<f> is 
convex and $ restricted to such set is an isometry. 

2.3. On the isometric dynamics. Whenever we have conjugate dynamics (i.e. a commuting 
diagram), we are allowed to carry certain properties from the classical to the quantum setting. For 
instance, we note that for A, and T given in the introduction, all being linear, are such that convex 
sets are taken to convex sets by these maps, and the same is true for its respective inverse images. 
Also supposing the commuting diagram it is straightforward to show that if / is a fixed point for T 
then Kfj,(f) is a fixed point for A. We have seen in propositions 12. 11 1 2.41 that in certain regions we have 
an isometric dynamics for A, but such fact is not a consequence of the classical behavior. Instead, 
this is intrinsic to CP dynamics, which are always non-expansive with respect to the trace distance. 

Example. Consider the 2-qubit CNOT operator, for which an asymptotic analysis is made in 
[13]. Consider $(p) = C\pC\ + CipCi, where Ci\i,j) = \i,i(Bj), C^l^j) = It can be shown 

that 



(2.6) 




m=0 



a 



c 



c 



c 



a 



c 



c 



c 



(2.7) 




c b d d 
c d b d 
c d d b 




c b d d 
c d b d 
c d d b 
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where p(k) := $ fc (p) and a,b,c,d are numbers depending on the initial state p only. A routine 
calculation shows that [a\ , cr 2 ] =0 and that the same holds for any two elements on the segment 
connecting G\ and a 2 (such segment is, of course, co{o~\, cr 2 }). 



More generally, does commutativity on co O(p) holds for every CP map considered, or there exists 
a dependence on the classical dynamics and/or the quantization choice? Motivated by the above 
example, we have that the following holds: 

Proposition 2.5. a) Suppose that the diagram U.S\) commutes for a certain choice of p and f E L°° . 
Ifg E L°° and f E Q(g) then Kfl (f) E tt A {^{g)). b) \p,rj\ = for every p,r, E ^0(0, VC E D(C N ). 



Proof, a) Immediate from the commutative diagram, b) Let p,r] E co 0(C). Then 

(2.8) p = limp m , p m Eco 0(C), p m = Y] X m ,iPm,i, Pm,i = um & m ' iJ (0 

m ' ' j 

i 

(2.9) r] = \imr] n , r] n EcoQ((), rj n = }X n ,krin,k, Vn,k = hm (() 

n I 

k 

Therefore, 

(2.10) p = lim V s A mii lim$ im ^ J (C), rj = lim A„ )fc lim $'™ fc - i (() 



n 

k 



The conclusion follows. 

□ 

Remark. It is not true that every two densities in the attractor set of a CP $ commute, as 
simple examples can show. However, if they belong to the convex hull of the same density matrix, 
then this is true as the above proposition shows. 



3. Quantization procedure 

3.1. Measures and a representation for M n . Measure theory has many applications to algebras 
appearing in quantum mechanics. Some measure-theoretic analysis that occur in quantum informa- 
tion problems are related to identifying properties of states. For instance, entanglement of formation 
can be described in terms of integrals over the state space [TTJ . In order to introduce our discussion, 
we recall the following definition. A barycenter of a probability measure p E Mi(K) over a com- 
pact, convex set K is an element b(p) E K such that tp(b(p)) = J dp(u), for every if) linear 
functional on K. Two basic results on barycenters of measures are the following [3]: 

Theorem 3.1. Let K be compact and convex, p E M\(K), the set of probability measures over K . 

a) There exists a unique point b(p) E K such that 

(3.1) </>(%)) = J ^Pdp, if)EA(K), 

where A(K) denotes the set of real continuous affine functions. 

b) There exists a net p a E M\(K) of measures with finite support which is weakly* convergent to 
p and is such that b(p a ) = b(p), for all a. 

Theorem 3.2. Let U be a C* -algebra and Eu the associated set of states. Then the set M^(Eu) of 
probability measures over E u with barycenter uj is a convex, compact subset of Mi(E u ), in the weak* 
topology. 
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We consider the C*-algebra M = M n of order n complex matrices. Recall that a state is u G M*, 
u : M n — > C such that it is positive (i.e., takes positive elements to positive real numbers, u(AA*) > 
0) and = 1. Let p G M n be a faithful state and identify its GNS triplet (H p , n p , Q p ) with 

(3.2) (C n g>C\M„® / n ,|VP» 
That is, 7T p (M n ) = M n <8> J n , where if p = J2i r i\ r i)( r i\ e then 

(3.3) \Vp) = J2 V*) ® I 7 "*) e C n ® C n 

i 

Also the associated state functional is Co> p : M n — >■ C, ^ P (X) = tr(Xp). It also holds that the linear 
span of vectors of the form \X^/p), X G M n , is dense in C <g) C And note that n p (M n )' — I n <g) M n . 

Remark. In principle we could have chosen a representation that takes L°° to M n , instead of 
M n ® /„,. We prefer to use the latter because this representation makes more explicit the relation 
between the finite algebra we are considering and its commutant. In a general setting, for a given 
C*-algebra A, k p takes L°° to the commutant of ir(A). This particular representation also appears 
in the modular theory for finite-level systems [2] . 

Example. (Orthogonal measures [3j [TTJ [17]) Denote by Em the state space for M = M n , the 
C*-algebra of matrices. Recall that Em is a convex and compact in the weak* topology. Denote by 
E{Em) the set of extreme points of Em- For p 6 Em, denote by M p (Em) the space of probability 
measures with barycenter p. Let u\ and U2 be positive linear functionals over M n and let u = uji+uj2- 
We say that u\ is orthogonal to u 2 , written u\ _L U2 if there exists a projection P 6 / n <E> M n such 
that 

(3.4) uj 1 {B) = (Pn w ,{B®I)n w ) 1 u 2 (B) = ((I-P)£l u ,(B®I)n u ) 

Let /ibe a positive regular Borel measure on Em- We say p is an orthogonal measure if 



(3.5) / uj dp{oj) _L / u dp{oj) 

Js Je m \s 

for every S C i?M Borel set. We note that orthogonal measures are in one-to-one correspondence 
with the abelian von Neumann subalgebras of the commutant tt p (M)' [3]. 

3.2. Transfer operator. With the notation of the previous section, let Em denote the convex 
weakly* compact set of states over M. Let M u (Em) be the set of probability measures over Em with 
barycenter uj. For p G M w (Em) define the map / G L°°(Em', p) n p (f) G ir p (M)', 

(3.6) (^,KJf)7rJA)n w ) = [ f(u')uj'(A) dp(u'), 



which is positive and non-expansive. We know by Tomita's theorem (3] that if p, is a nonnegative 
regular orthogonal Borel measure on Em then the map / t— > n p (f) is a ^-isomorphism of L°°(Em', p) 
into 7T P (M)' and 

(3.7) A = {K p {f):feL°°(E M - 1 p)} 

is an abelian von Neumann subalgebra of 7r p (M)'. Now define L : L°°(Em] p) — > L°°(Em] p), 
(3-8) L(f)(A) = J2^m(A)), 

i 

where F,i : M n — > M n , A G M n , f G Em, Pi > 0. Such operator is known as Ruelle trans- 
fer operator, or Ruelle-Perron-Frobenius operator. The spectral theory of this map is known as 
thermodynamic formalism [TJ. 



Remark. When studying transfer operators, it is usual to consider them acting on the space 
C(Em) of continuous functions. The analysis made in this article also works in that setting, except 
for the statements related to orthogonal measures and algebra commutants. Our preference for the 
more general space L°° is due to the fact that in such space one benefits from certain algebraic 
properties encountered when studying the operator k p in full generality, see [3J. 

While the interplay between the transfer operator L (acting on continuous functions) and its dual 
L* (acting on measures) has been extensively studied PQ, little is known about the behavior of the 
noncommutative counterparts of the measures (the state operators) with respect to quantum ana- 
logues of the transfer operator (acting on matrix algebras). We aim to describe a relation between L 
acting on L°° functions and completely positive maps (quantum channels) acting on a matrix alge- 
bra. Note that L itself is completely positive, since it is a positive operator acting on a commutative 
C*-algebra. 

3.3. Conditions on the CP map. We are interested in the evolution of L(f) and its relation to 

L°°(Em; /u) — ^ L°°(E M ;ti) — ^ ••• 

I n ® M n — I n <g> M n — • • • 
where A is a CP map to be determined. 

Definition. We say that k p is a CP quantization for the transfer map T : L°° — > L°° if 
there exists A CP such that the above diagram commutes. We say the quantization is faithful if 
/ G i — y K>n(f) G In ® M n is faithful as a representation. It is well-known that a quantization 
is faithful if and only if fi G S(M uj (Em)) [3]. 

We are interested in finding conditions over A in order to obtain a commutative diagram, that is, 
Kfj, o L = A o Ky,. We will suppose that A : I n <g> M n — > I n <g> M n is completely positive. Note that 
I®M n ~I n ® B(C n ) C B(C n ® C n ) ~ B(C n2 ). We will then take VJ : C n ^ C n such that 

(3.10) A(/ n g) x) = ^(/ n g) ® x)(/ n <8» y;) 

i 

We have 

(3.11) K„(L/)7r p (A)n p ) = (/ ® M Lf )ix p {A)n p ) = ! J2Pif(FMW(A) dfi(u'), 

J i 

where above we write n^Lf) — I (g) M L f. On the other hand, 

(Q p ,A(^(f))n p (A)Q p ) = ({l p ,A(I®M f )n p (A){l p ) 

(3.12) = (n p ,^(I n <g>VMI n ®M f ){I n ®V;)Tr p {A)n p ) = (fi p ,^/ n ®(V i M / V7')7r p (A)n p ) 

i i 

Therefore we are interested in conditions over the in order to obtain 

(3.13) I n ®M Lf = I n ®Y< V * M f V * 

i 

In order to find a solution to the above equation, it is enough to solve 

(3.14) M Lf = Y i V i M f V; 

i 
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(3.9) 



Now we perform the general procedure and for clarity the reader may consider at each step the 
case n = 2. Let \x G M p (Em„), the set of measures over the state space for M n , with barycenter p. 
Let / G L°°(/i). Then K^(f) G 7r p (M n (C))' = /„, Cg> M n (C). It is a simple calculation to show that 
{7r p (£'j :; )|y / p)}j J= i i .... n is a basis for C n Cg> C n , where are the matrix units for M n . For instance, if 
n = 2, write |0) = (1,0), |1) = (0,1), \^p) = p \0) ® |0) ® |1). Then 

TTp^n)!^} = (^ii <S) J 2 )(p |0> <S) |0> +pi|l> O |1» 

(3.15) = po(£n ® I 2 )(|0) ® |0)) + V i{En ® h){\l) ® |1» = (po, 0, 0, 0) T 
Similarly, 

(3.16) 7r p ( J E; 12 )|Vp) = (0,p 1 ,0,0) T , 7r p ( J E 21 )|Vp) = (0,0,p ,0) T , tt p (£ 22 )| y'p) = (0, 0, 0, Pl ) T 

To avoid too many indices we make a pair correspond to the number n(i — 1) + j. That 
is 11 — )- 1, 12 — 2, 21 — > 3, 22 — > 4 and so on. In this way we relabel the matrix units. Denote 
Vk = ^ p{Ek)\\fp) ■, for k = 1, ... , n 2 . In such basis we have 

(3.17) [«„(/)]„ = (viMffo) = (Vpk P (^)ls(/)^(^)IVp) 

(3.18) = (Vp|^(/)7r p (^)7r p (^.)iVp) = (VpK(/K(^)IVp) = / f^'WiElE^iu') 



A routine calculation proves that the matrix K^(f) is block diagonal and all blocks are equal. For 
instance, if n = 2 such matrix has order 4 and 



TP* TP TP TP* TP 



E* 2 E 2 



E 2 , E{E 3 
E 4 



0, E* X E A = 



E^E 3 = 0, i? 2 i?4 = 



E*E 1 = 0, E* 3 E 2 = 0, E*E 3 
E\E\ = 0, E\E 2 = 0, EIE3 



Ei, E3E4 
E 3 , E\E± 



E 2 
E 4 



Hence we obtain the matrix 

7 / uj{E x )dn J f u(E 2 )dfi 

j f u(E 3 )dfi ffu(E A )dfi 

/ / u{Ei)dfj, j f oj(E 2 )dfi 

jfu(E 3 )dfi ffu(E 4 )dfi 



(3.19) 



Jfu{E x )diJL ffu(E 2 )dn 
Jfoo(E 3 )dp Jfoo(E 4 )dfi 



I <g> M 



thus providing an explicit expression for M/, and analogously for Mr,f. So equation (13. 14ft becomes, 
for n = 2, 



M Lf = J2ViM f V* & 



(3.20) Y,P> 



J f(F^ uiE^dfi J f(Fi) u(E 2 )d/i 
J f(F t ) u(E 3 )dfj, J f(Fi) u(E A )dfj, 



i=l 



'Jfu}{E x )dn Jfoo(E 2 )dp 
jfu{E 3 )dfi Jfu(E 4 )dfx 



V* 



The number is determined by the transfer operator and can be at most n 2 , a restriction obtained 
from the Kraus decomposition. To solve the above equation, it is enough to look for V, which are 
diagonal, and we also set I = k. For a general matrix A = ((%), we have 

(3.21) Vi = diag« ...,<) V^F/ = fej, r, s = 1, . . . , n, i = l,...,k 



Now assume a correspondence between the z-th matrix on each side of equation (13.201) . Then for the 
diagonal Vi above and 

(3.22) A = M f = [J fu(E ra )dn], r,s = l,...,n 
we get from (I3.20p . 

(3.23) pi[ J f(Fi) ui(E rs )dp] = v l r vi J f u)(E ra )dfi, r, s = 1, . . . , n, i = 1, . . . , k 
For n = 2, this becomes 



J f{Fi) uiEJdfi, j f{Fi) u(E 2 )dfi 
J f(E) ou(E 3 )dfi J f(Fi) u(E 4 )dv 



(3.24) Vl 

Note that the equations obtained from the diagonal are 



ii 2 //^(£i)^ 

ffu(E 3 )d» 



\v 



v{v l 2 f f u)(E 2 )d(i 



A 12 



/ / u(E 4 )dn 



(3.25) 
and so 

(3.26) 



Pi / f(Fi(uj))uj(E rr )dfi(u) = \v l r \ 2 / f(u)u)(E rr )dp,(u>), i = 1, . . . , k 



I Pi J f{F i {u))u{E rr )d l i{u) 
j f(u)u{E rr )dfi{uj) 



,n, 



k 



On the other hand the equations obtained from non-diagonal entries must also satisfy such conditions, 
but this is not expected to hold in general. So we have restrictions on the states with respect to a 
given measure \i. 

Example. Take n = 2. If we choose \i — 5 Uo , we obtain the equations 

Pi f(Fi{u )) = \v[\ 2 f(u ) 



(3.27) 
(3.28) 
(3.29) 

So if f(u)o) ^ we conclude 
(3.30) 



pJ(F t (co ))=v\vif(u ) 
p t f(E(u )) = K| 2 /M 



Kl 2 



f(EM) . 



k 



Note that if ojq is a fixed point for the CP map A : p i— > ^2iPiFi then 
(3.31) v i = vi = ^/Pie ie % i = 1,2, f e R 

This establishes the existence of a solution for the commutative diagram described above in the 
present case. Note that if YliiPi = 1> as ^ i s usually assumed, then A(J) = V{V* + V 2 V 2 * = I so we 
have that if ojq is a fixed point for A which is mixed-unitary then each Fi is fixed by such point and 
therefore A is unital and trace-preserving (see section [2TTj) . 

An open question is to find all solutions for a general quantum channel, i.e., in the case that not 
all Vi are diagonal. Also the case in which /x is a convex combination of Dirac deltas is of interest, 
a problem where a generalization of the above method seems possible. Another open question is to 
describe the solution when we consider an orthogonal measure. 
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Example. Consider M 2 (C), let p, = XiS^ + \2Su2, for u\, U2 fixed states with support on the 
diagonal units of M2. Write n^Lf) = I ® Ml/, K/j,(f) — I ® Mf. Then we wish to obtain V^- 
satisfying = VjMfV*, and again we seek V} which are diagonal. For each i, we have: 

(3.32) pi J f(Fi(u))) u{E rs )dp = v l T v\ J f(u) u(E rs )dp, r, s = 1, . . . , n, i = l,...,k 

We conclude that 

p i (x 1 f(F i (u l ))u l (E n )+X 2 f(F i (u 2 ))oj 2 (E ll ] 



(3.33) 



(3.34) 



Kl 2 



\v. 



i 12 



XlfiiO^LO^Eu) + X 2 f(w2)u2(Eu) 
P AX l f{F l {u l ))u l {E 22 ) + X 2 f{F i {uj 2 ))u J2 {E22] 



2 Xif(u 1 )u 1 (E 22 ) + A 2 /(w 2 )w2(-E , 22) 

The above reasoning can be generalized and summarized in the following: 



Theorem 3.3. Let L : C(E Mn ) -> C(#mJ, L(g)(p) = Ei=iPi9(Fi(p)) , k < n 2 , Pi > 0. 

(1) For a given f > and for every p = J2 m ^mi m , X m > 0, Ui with support on the diagonal 
units of M n (C), there exists a solution for A o = K^Lf), given by A(p) = Xa=i ViPVi > 

Vi = diag(v{, . . . ,v l n ), and 



(3.35) 



Pi{ Em X m f(Fi{uJ m ))u m {E. 



J J 



\ Y^ m ^mf{Um)Um{Ejj] 



l,...,k, j = l,...,n 



(2) If Lf = J2iPif(Fi) with Fi = Ui ■ U* , Ui unitary, there exists Vi which satisfies A o = 
Kfj,(Lf) for every f G L°° if and only if the Ui are all fixed points for $(p) = YliPi-^i(p) ■ 

4. Classical diagram and open questions 

We recall that a density matrix p is a fixed point for a mixed unitary channel A if and only if 
such matrix is the barycenter of a measure which is invariant for the associated Markov operator P. 
If A(p) = ^ZiPiFi(p), Fi(p) = UipU* then we know that A(p ) = Po implies that Fi(p ) = p for each 
i [13]. The associated Markov operator Pm on measures is given by 



(4.1) 



P M p{B) =Y] Pi dp 



Classical Markov operators are related to the discussion in this paper in the following way. Consider 
the following diagram: 

Pm , 



(4.2) 



L,P 



L°°(E M ;n) — + L°°(E M ;p) 



The lower part of the diagram is the one considered in the previous sections. Above, M a denotes 
the family of absolutely continuous (AC) measures, / is the integral formula map / i->- pf, Pf(A) = 
j A f dp and R stands for the Radon- Nikodym derivative of an AC measure. Assume that Pm is a 

Markov operator on measures acting on absolutely continuous measures. Now let / G L 1 , f > and 
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consider p,f. Then Pkf(p/) is also AC, so we can write it in the form PMPf{A) = f A g(p) dfi(p), where 
g is the Radon- Nikodym derivative with respect to the Borel measure Pm^j on D(C n ). Therefore 
for every / G L 1 we obtain a unique g G L 1 , such correspondence is a Markov operator P acting on 
densities (positive functions) and it is such that the upper part of the above diagram (not involving 
matrices) commutes [7]. 

On the other hand, suppose we are given A which can be written as A(p) = J2iPi^i(p)- Then we 
can obtain a Markov operator Pm on measures, via equation (14. ip . and therefore another operator P 
acting on densities. It is an open question to determine in which cases we have that the P obtained 
is such that the lower diagram commutes. This is the inverse problem: given a CP map on a matrix 
algebra, obtain a classical approximation. 

Another open question is the following: given L operator on positive functions and / a function, 
consider a quantization A f Markov operator Pm on measures and by the above construction 

an operator P on positive functions. Then, we wish to determine how similar are L and P, and 
whether we can find cases in which L = P (or at least Lf = Pf), and whether such solutions are 
nontrivial. 

5. Appendix 

In this appendix, we present the results on non-expansive maps cited in this work. The main ideas 
have been presented in [UEHH], and are adapted here for our purposes. For completeness we present 
the proofs, which are mostly slight variations of the works mentioned. 

Proposition 12.11 [I] For all p G D(C ) we have the following, a) f2(p) is minimal under $. b) 
For every k > 1, $ fc is an isometry on fi(p). c) fl(p) is strongly invariant under $. 

Proof, a) Let 77 G f2(p), rj = linn^oo $ rafe (p), n k — > oo as k — > oo. For any I > 1, $'(77) = 
lim^oo <5> nk +\p) g Q(p). Hence 7(77) C Q{p) = Q{p). Now let ( G Q(p), C = hm*^ $ Tfc (p), n -»■ 00 
as k — > 00. Without loss of generality assume s& := t& — > A;, k — 1, 2, . . . . Now note that 

(5.1) d($'*faU) < d(<S> s *( V ),$ s i° +n *(p)) + d($ T *(p),() < d( V ,$ n *(p)) + d($ T *(p),() 

so we deduce that ( G fi(p). Hence we have fl(p) C 0(77) C 7(77) and so we get = fi(?7) = 7(77). 
Therefore fi(p) * s minimal. 

b) Note that from a), we have that fl(p) = Q{r]) so there is {<7k} C N, o"fc — >■ 00 such that 77 = 
limbec <& ak (rj). We claim that $ CTfe (C) -» £, A; -» 00, for any ( G O(p). Indeed, if C = limfc->oo $ Sft (?7) 
then 

d($ CTfc (C),C) < d($ fffc (C),^ fc+Sfc (^)) +d(<S> ak+Sk (7]),<S> Sk (ri)) +d(<S> Sk (r])X) < 

(5.2) < 2d(C,^ Sk (v)) + d(^ ak (v),v) -> 

as A; — >■ 00. In the last inequality we have applied the non-expansiveness of $ and also due to this 
fact we have for any k > 1, 

(5.3) ^(77,0 = lim d($^(?7),<£ afc (C)) < d($ fe (77),$ fe (C)) < d(77,C) 

k—>oo 

This shows that $ fc is an isometry on fi(p). 

c) From b) we get that $ fc : f2(p) — > &(p) is one-to-one. To prove that $ fc is onto f2(p), we fix 
C G fi(p). For large /, m, n, 

d($ an - k ((),<$> am - k (()) < 
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(5.4) < 2d($ ,Ji (C), C) + d($ ff " +a ^ fc (C), ® am+ai - k (()) 

Also from b) we have that $ CTi (C) — ^ C ^ — ► °°- Moreover, for fixed 1, {$ an+a i~ k (()} is Cauchy, 
converging to $ <Ji_fc (C) as n — > oo. It follows that {$ cr ' l_fc (C)} is Cauchy, so let 77 denote its limit. 
Then it is clear that $ fc (ry) = lim^oo 3> <Tn (£) = C, so <3> fc is onto. This proves the strong invariance of 

n(p). 

□ 

Lemma 12.21 [5] If the restriction of $ to a subset i? of D(C ) is an isometry then the restriction 
of $ to co B is also an isometry. 

Proof. To prove that d($(p), $(77)) = d(p, 77), for all p, 77 G co B, it is enough to prove the equality 
for pi E co B and p G B: 

(5.5) d{${p),$(ri)) = d{p,rj), p 1 ecoB,r)EB 

In fact if this is true then <3> restricted to Sx = B U {pi} is an isometry so the above equality holds 
for pi G B\ (but also p\ G co B) and p 2 G co B\ = co(B U {pi}) = co £>. To prove (15.51) we use 
induction on the positive integer k in the representation p = Yli=i \Pii \ > 0, Yli=i — 1> Pi ^ 
For fc = 1, the claim is true by assumption. Suppose it is true for k — 1 and write 

fc-i 

(5.6) p = ^ XiPi + A fc p fe = api + (1 - a)p 2 

i=i 

where 

fc-i ^ fc-i 

(5.7) a = y]Ai, pi = -Y^AiPi, p 2 = Pfc 

^ a ^ 

4 = 1 1=1 

Now note that p 2 G 5 and pi satisfies (15.51) by induction hypothesis. Therefore the triangles f]pip 2 
and $(r7)$(p 1 )$(p 2 ) are congruent. In particular <i(<3>(pi), $(p 2 )) = ^(pi,p 2 )- Since p belongs to the 
segment generated by pi and p 2 we have 

(5.8) d{p, Pl ) = d{p h p 2 ) - rf(p, p 2 ) < rf($(pi), $(p 2 )) - $(p 2 )) < rf($ (pi), $(p)) < p) 

It follows that rypip and $(?7)$(pi)$(p) are congruent. So (15. 5p holds and this concludes the lemma. 

□ 

Proposition [2731 [4 J a) The set co fi(p) is strongly invariant under $. Also there is an affine group 
of isometries on the closed linear variety spanned by fi(p) which coincides with $ on coQ(p). b) 
co fl(p) contains exactly one fixed point p of $, given by 

k 

(5.9) p = lim ~ V$ m (?7), VpG^oT^) 

fc^oo /c z — ' 
m=0 

Proof, a) The first claim follows from a convexity reasoning over fi(p), which is strongly invariant: 
the surjective argument is immediate and injectivity follows from the fact that $ fe is an isometry 
on co fi(p) also, by lemma 12721 Let L denote the affine space spanned by fi(p). For po G L, write 
p Q = At7+(1-A)C, rp C e co n{p). Define T k {p ) := \<5> k (r]) + (l- A)$ fe (C), fc > 0. It is straightforward 
to check that T k (p ) does not depend on the particular convex expression for p . Also it is clearly 

affine, satisfying T° — I and T p+q = T p T q , p,q>l. By continuity we extend the definition to L. 
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b) Suppose 770 G co Q(p) is a fixed point for $. Let £ G and let £ G co 7(£), say £ 

EILi Ai$ fli (0, < Ai < 1, Ei = 1, Si > 1, i = 1, . . . ,n. We have 

m=0 i=l m=Si 

= E H [ E $m ^) - E *'(o + E $m (0 

i=l m=0 i=0 m=fc 

(5.10) ~ $m (o + E A 4 [ - E *'(o + E $m (^) 

m=0 i=l i=0 m=k 

Hence since 7(£) is bounded, 

T]0 



hmsupd(i ® m (0,Vo) = limsupd(± ]T $ m (C),' 



(5.11) =limsu P( i(^$ m (C),$ m (r/o))<rf(C,%) 



m=0 



But fi(p) is minimal so we get that ( G co 7(£) = co fi(p) and 770 G co fi(p). Therefore 

(5.12) inf d(C,Vo) = 

CGco 7 (0 

So we have r/o = lim^oo t Em=o ^* m (0 f° r ever y £ £ ^(p) an d therefore for every £ G co fi(p). In 
particular r/ is the unique fixed point for $ in co fi(p). 

□ 

Proposition l2~T4l a) There exists a non-expansive projection r onto b) We have that Q$ is 
convex and $ restricted to such set is an isometry. 

The proof uses the following lemma, which can be seen in [8]. 

Lemma 5.1. Every $ : D(C N ) — > D(C N ) CP map is such that every subsequence of {$ fc } has a 
convergent subsequence that converges uniformly on compact subsets of D(<C N ). 

Proof of proposition 12.41 a) By the above lemma there exists a subsequence {$ fci } of {$ fc } that 
converges uniformly on compact subsets of D(C N ). By passing to a subsequence and applying the 
lemma —7- 00 as j — » 00, and {$ ^} converges uniformly 

on compact subsets of D(C N ). Define r : D(C N ) D(C N ), 

(5.13) r(p) := lim $^ +1 "^(p), p G D(C N ) 

To see that r is a projection onto first note that r(p) G f2$ for every p G D(C M ) so we have 
t(D(C n )) C Let 77 G It follows from proposition 12.11 that 77 G Q(r]) and as 7(77) has a 

compact closure, $(fi(r/)) = Therefore there exists for each k > 1 a point 77^ G ^(77) such that 

Q k {Vk) — V- Denote by g the pointwise limit of {$ fc,; } and let {77^} be a convergent subsequence of 
{Vki} with limit ( G Q(rj) (such a subsequence exists as fl(r)) is compact). Then note that 

(5-14) d(^( Vlj ),g(C)) < d(^( Vh ),g( Vlj ))+d(g( Vlj ),g(0) -+ 0, 
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as j — > oo, since } converges uniformly to g on fl(r]). Therefore g(() = r]. Now note that 

(5.15) = ($^ (£)) 

for each i > 1. The LHS converges to g(() whereas the RHS converges to r(g(()). Therefore r(pj = r] 
for all r\ G so r is a non-expansive projection onto 

b) We have for all p,r/ <E 

(5.16) d(p,r?) = d(T(p),r(ri)) < lim d($^ +1_ ^(p), ^ ki i+ 1 ~ ki i ( v )) < d($(p), $(?/)) < d(p,r/) 

Let r : D(C N ) — >■ f2$ be the projection given in a). For every p, r? G f2$ and < A < 1 we have 

(5.17) Ap + (1 - A)r7 = Ar(p) + (1 - A)r(^) = r(Ap + (1 - A)??) G 

□ 
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